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The part of the proton spin E carried by u, d, s quarks is calculated in the framework of the QCD 
sum rules in the external fields. The operators up to dimension 9 are accounted. An important 
contribution comes from the operator of dimension 3, which in the limit of massless u, d, s quarks 
is equal to the derivative of QCD topological susceptibility x'(0)- The comparison with the experi- 
mental data on E gives x'(0) = (2.3 ± 0.6) x 10~ 3 GeV 2 . The limits on E and x'(0) are found from 
selfconsistency of the sum rule, E > 0.05, x'(0) > 1.6 x 10~ 3 GeV 2 . The values of g A = 1.37 ± 0.10 
and g\ = 0.65 ± 0.15 are also determined. 
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' In the last years, the problem of nucleoli spin content and particularly the question which part of the nucleon 
t-H , spin is carried by quarks, attracts a strong interest. The valuable information comes from the measurements of the 
spin-dependent nucleon structure functions g\ (x, Q 2 ) in deep inelastic e(fi)N scattering (for the recent data see [1,2], 
for a review ||). The parts of the nucleon spin carried by u, d and s-quarks are determined from the measurements 

■ of the first moment of gx(x, Q 2 ) 

T P ,n(Q 2 ) = J dxg 1;p , n (x,Q 2 ) (1) 


The data allows one to find the value of E - the part of nucleon spin carried by three flavours of light quarks 
. E = Au + Ad + As, where Am, Ad, As are the parts of nucleon spin carried by u, d, s quarks. On the basis of the 
y—{ ' operator product expansion (OPE) E is related to the proton matrix element of the flavour singlet axial current j® 5 

o : 

«?■ 2ms^=(p,s\jl 5 \p,s), (2) 

On ; 

where is the proton spin 4-vector, m is the proton mass. The renormalization scheme in the calculation of 
fS { perturbative QCD corrections to r Pi „ can be arranged in such a way that E is scale independent. 

An attempt to calculate E using QCD sum rules in external fields was done in rcf. Let us shortly recall the 
idea. The polarization operator 

43 : 



X 



n(p) = i J d?xe vx (0\T{r](x),fj(0)}\0) (3) 

was considered, where 

Vfr) = c abc U°(*)<?7X(^ 1^5d c (x) (4) 

is the current with proton quantum numbers [5], u a , d b are quark fields, a, b, c are colour indeces. It is assumed that 
the term 

AL = J° 5 ^ (5) 

where is a constant singlet axial field, is added to QCD Lagrangian. In the weak axial field approximation II (p) 
has the form 

U(p)=U^(p) + U^(p)A,. (6) 

nji (p) is calculated in QCD by OPE at p 2 < 0, \p 2 \ R~ 2 , where R c is the confinement radius. On the other 
hand, using dispersion relation, II^ (p) is represented by the contribution of the physical states, the lowest of which 
is the proton state. The contribution of excited states is approximated as a continuum and suppressed by the Borcl 
transformation. The desired answer is obtained by equalling of these two representations. This procedure can be 
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applied to any Lorenz structure of H^' (p) , but as was argued in [6,7] , the best accuracy can be obtained by considering 
the chirality conserving structure 2p fJ _pj 5 . 

An essential ingredient of the method is the appearance of induced by the external field vacuum expectation values 
(v.e.v). The most important of them in the problem at hand is 

(o|j£ 5 |o) = 3/ % (7) 

of dimension 3. The constant /q is related to QCD topological susceptibility. Using (5), we can write 

<0|j« 5 |0).4 = lim q ^ i J d 4 x e ^(0|Ta 5 ( s ),^ 5 (0)}|0}^ = 

= lim q -,oP flu (q)A v (8) 

The general structure of P^ v (q) is 

P„u(q) = -P L (q 2 )S^ + P T (q 2 )(-S^q 2 + q^q v ) (9) 

Because of anomaly there are no massless states in the spectrum of the singlet polarization operator P M „ even for 
massless quarks. Pr,L(q 2 ) also have no kinematical singularities at q 2 = . Therefore, the nonvanishing value P^(0) 
comes entirely from P^(g 2 ). Multiplying P fll/ (q) by q^q^, in the limit of massless u,d,s quarks we get 

q„quP„u(q) = -P L {q 2 )q 2 = N 2 (a s / Avfi J d^xe^x 

x (0\TG^(x)G^(x),GZ(0)GZ(0)\0), (10) 

where G™„ is the gluonic field strength, G M „ = (1 /2)e ^vXaGxa .(The anomaly condition was used, Nj = 3.). Going to 
the limit q 2 — > 0, we have 

f 2 = -(l/3)Pfc(0) - ^^'(0), (11) 
where x(<? 2 ) is the topological susceptibility 

X {q 2 ) = i /d 4 .T e ^(0|Tg 5 ^),Q 5 (0)|0) (12) 



Q 5 (x) = K/8tt) G^(x)G^(0), (13) 

As is well known (see, e.g., the review S), x(0) = if there is at least one massless quark. The attempt to find x'(0) 
itself by QCD sum rules failed: it was found Q that OPE does not converge in the domain of characteristic scales for 
this problem. However, it was possible to derive the sum rule, expressing E in terms of f§ (7) or x'(0). The OPE up 
to dimension d = 7 was performed in ref. [|J. Among the induced by the external field v.e.v. 's besides (7), the v.e.v. 
of the dimension 5 operator 

.9(0|^g7a(l/2)A"GS /3 g|0) A = 3/ l0 A /3 , q = u,d,s (14) 
i 

was accounted and the constant ho was estimated using a special sum rule, 

ho w 3 x 10~ 4 GeV i . There were also accounted the gluonic condensate d = 4 and the square of quark condensate 
d = 6 (both times the external field operator, d — 1). However, the accuracy of the calculation was not good 
enough for reliable calculation of E in terms of /q : the necessary requirement of the method - the weak dependence 
of the result on the Borel parameter was not well satisfied. 

In this paper we improve the accuracy of the calculation by going to higher order terms in OPE up to dimension 
9 operators. Under the assumption of factorization - the saturation of the product of four-quark operators by the 
contribution of an intermediate vacuum state - the dimension 8 v.e.v. 's are accounted (times A^): 

- g(O\qa a0 (l/2)X n G^q • qq\0) = m 2 (0\qq\0) 2 , (15) 
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where TOq = 0.8 ± 0.2 GeV 2 was determined in S. In the framework of the same factorization hypothesis the induced 
by the external field v.e.v. of dimension 9 



a s (0\j^\0)A{0\qq\0y 



(16) 



is also accounted. In the calculation we used the following expression for the quark Green function in the constant 
external axial field H: 



(0|T{«£(a;), $(0)}|0) A = iS ab x a p/2ir 2 x 4 



+ (l/27T 2 )S ab (Ax)( l5 x) aP /x i - (l/12)S ab S aP (0\qq\Q) + 



+(l/72)iS ab (0\qq\0){xA l5 - Ax l5 ) al3 + 



(1/12) f£5 ab (A~f 5 ) a p + (l/216)S ab h 



(5/2)x 2 A l5 - {Ax)x~f 5 



(17) 



a/3 



The terms of the third power in x-expansion of quark propagator proportional to A^ are omitted in (17), because 
they do not contribute to the tensor structure of II M of interest. Quarks are considered to be in the constant external 
gluonic field and quark and gluon QCD equations of motion are exploited (the related formulae are given in [10]). 
There is also an another source of v.e.v. ho to appear besides the x-expansion of quark propagator given in eq.(17): 
the quarks in the condensate absorb the soft gluonic field emitted by other quark. A similar situation takes place also 
in the calculation of the v.e.v. (16) contribution. The accounted diagrams with dimension 9 operators have no loop 
integrations. There are others v.e.v. of dimensions d < 9 particularly containing gluonic fields. All of them, however, 
correspond to at least one loop integration and are suppressed by the numerical factor (2ir)~ 2 . For this reason they 
are disregarded. 

The sum rule for £ is given by 



E + C M 2 = -1 



9A 2 

N 



2 /M : 



{a 2 L 4 / 9 + 



+6n 2 f 2 M* El (^-\ L-^ + Un 2 h M 2 E 1 1 




(18) 



Here M 2 is the Borel parameter, Ajv is defined as = 32tt 4 X 2 v = 2.1 GeV 6 , (0\ri\p) = Xnv p , where v p is proton 
spinor, W 2 is the continuum threshold, W 2 = 2.5 GeV 2 , 



a = -(2tt) 2 (0|to|0) = 0.55 GeV 3 



(19) 



E (x) = l-e~ x , E^x) = 1 - (1 + x)e~ x 

L = ln(M / A)//n(/i/ A), A = Aqcd — 200 MeV and the normalization point u was chosen /i = 1 GeV. When 
deriving (18) the sum rule for the nucleon mass was exploited what results in appearance of the first term, -1, in 
the right hand side (rhs) of (18). This term absorbs the contributions of the bare loop, gluonic condensate as well 
as a s corrections to them and essential part of terms, proportional to a 2 and ra^a 2 . The values of the parameters, 
a,}? N ,W 2 taken above were chosen by the best fit of the sum rules for the nucleon mass (see [11], Appendix B) 
performed at A — 200 MeV . It can be shown, using the value of the ratio 2m 8 /(m u + m<i) = 24.4 ± 1.5 [12] 
that a(l GeV) = 0.55 GeV 3 corresponds to m s (l GeV) — 153 MeV. a s corrections are accounted in the leading 
order (LO) what results in appearance of anomalous dimensions. Therefore A has the meaning of effective A in LO. 
The unknown constant Cq in the left-hand side (lhs) of (18) corresponds to the contribution of inelastic transitions 
p — > N* — > interaction withA^ — > p (and in inverse order). It cannot be determined theoretically and may be found 
from M 2 dependence of the rhs of (18) (for details see [11,13]). The necessary condition of the validity of the sum rule 
is |S| 3> \CoM 2 \exp[(— W 2 + m 2 )/M 2 ] at characteristic values of M 2 [13]. The contribution of the last term in the 
rhs of (18) is negligible. The sum rule (18) as well as the sum rule for the nucleon mass is reliable in the interval of 
the Borel parameter M 2 where the last term of OPE is small less than 10 — 15% of the total and the contribution of 
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continuum does not exceed 40 — 50% . This fixes the interval 0.85 < M 2 < 1.4 GeV 2 .The M 2 -dependence of the rhs 
of (18) at /q = 3 x 10 -2 GeV 2 is plotted in Fig.l. The complicated expression in rhs of (18) is indeed an almost linear 
function of M 2 in the given interval! This fact strongly supports the reliability of the approach. The best values of 
E = S /H and C = C /lt are found from the % 2 fitting procedure 

1 " 

X 2 = ' E ~ C * M i ~ i? ( M ')l 2 = ™ m < ( 2 °) 

»=i 

where R(M 2 ) is the rhs of (18). 

The values of E as a function of f§ are plotted in Fig. 2 together with vx 2 . In our approach the gluonic contribution 
cannot be separated and is included in E. The experimental value of E can be estimated [1,2] (for discussion see [|14[) 
as E = 0.3 ± 0.1. Then from Fig.2 we have / 2 = (2.8 ± 0.7) x 10~ 2 GeV 2 and x'(0) = (2.3 ± 0.6) x 10~ 3 GeV* . 
The error in f§ and x' besides the experimentall error includes the uncertainty in the sum rule estimated as equal 
to the contribution of the last term in OPE (two last terms in Eq.18) and a possible role of NLO a s corrections. 
At /g < 0.02 GeV 2 \ 2 is much worse and the fit becomes unstable. This allows us to claim (with some care, 
however,) that x'(0) > 1.6 x 10~ 3 GeV 2 and E > 0.05 from the requirement of self consistency of the sum rule. The 
X 2 curve also favours an upper limit for E < 0.6. At f 2 = 2.8 x 10~ 2 GeV 2 the value of the constant Co found 
from the fit is Co = 0.19 GeV~ 2 . Therefore, the mentioned above necessary condition of the sum rule validity 
is well satisfied. Recently, the first attempt to calculate x'(0) on the lattice was performed jn|. The result is 
x'(0) = (0.4 ± 0.2) x 10~ 3 GeV 2 , much below our value. However, as mentioned by the authors, the calculation has 
some drawbacks and the result is preliminary. 

From the same sum rule (18) it is possible to find g\ - the proton coupling constant with the octet axial current, 
which enters the QCD formula for r Pi „ 0. There are two differences in comparison with (18): 

1. Instead of f 2 it appears the square /| of the pseudoscalar meson coupling constant with the octet axial current. 
In the limit of strict SU(3) flavour symmetry it is equal to f 2 , f v — 133 MeV. However, it is known, that SU(3) 
symmetry is violated and the kaon decay constant, Jk ~ 1-25/^ . In the linear in s-quark mass m s approximation 
f v = 1.31/^. We put for /| the value /| = 2.6 x 10~ 2 GeV 2 , intermediate between f 2 and f 2 . 

2. ho should be substituted by m 2 f 2 . The constant m\ is determined by the sum rules suggested in Jl(|. A new fit 
corresponding to the values of the parameters used above, was performed and it was found; m 2 = 0.16 GeV 2 . 

The M 2 -dependence of g\ + C$M 2 is presented in Fig.l and the best fit according to the fitting procedure (20) at 
1.0 < M 2 < 1.3 GeV 2 gives 



5 ^ = 0.65 ±0.15, C 8 = 0.10Ge^ 2 y/x 2 = 1.2 x 10~ 3 (21) 

(The error includes the uncertainties in the sum rule as well as in the value of /|). The obtained value of g\ within 
the errors coincides with g\ = 0.59 ± 0.02 JlTj found from the data on baryon octet /3-decays under assumption of 
strict SU(3) flavour symmetry and contradicts the hypothesis of bad violation of SU(3) symmetry in baryon axial 
octet coupling constants Q . 

A similar sum rule with the account of dimension 9 operators can be derived also for g& - the nucleon axial /3-decay 
coupling constant. It is an extension of the sum rule found in [6] and has the form 



gA + C A M 2 ^l + ^ r e m2 / M2 
9\~ N 



2 r 4/9 i r, 2 2y2 7i,r2 ^ 2 m . ^ 2 a 

a L + 2tt mJ n M - -a— + o^UtTo 



4 M 2 3 JTr M 2 



(22) 



The main term in OPE of dimension 3 proportional to f 2 occasionally was cancelled. For this reason the higher order 
terms of OPE may be more important in the sum rule for g A than in the previous ones. The M 2 dependence of 
gA — 1 + C4M 2 is plotted in Fig.l, lower curve; the curve is almost the straight line, as it should be. The best fit 
gives 



g A = 1.37 ± 0.10, C A = -0.088 GeV' 2 , vV = 1-0 x 10~ 3 (23) 

in comparison with the world average gA = 1.260 ± 0.002 [19]. The inclusion of dimension 9 operator contribution 
essentially improves the result: without it gA would be about 1.5 and x 2 would be much worse. 
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FIG. 1. The M^dependence of £ + C Q M 2 at / c f = 3 x KT 2 GeV , eq.18, g\ + C 8 M 2 , and g A - 1 + C A M\ eq.22 
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